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Nuclear dynamics of giant resonances are investigated with the real-time Skyrme TDHF 
method. The TDHF equation is explicitly linearized with respect to variation of single- 
particle wave functions. The time evolution of transition densities are calculated for giant 
dipole resonances. The time-dependent densities of protons and neutrons suggest that 
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the dynamics of giant dipole resonance in neutron-rich nuclei are significantly different 
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from that in stable nuclei with N « Z. 



1. Introduction 



Atomic nuclei exhibit a variety of collective modes of excitation. In particular, the giant 
resonances always have been of central interest in nuclear structure and reaction studies. 
The giant resonances correspond to the most fundamental oscillations in nuclei, in the 
sense that they exhaust major portions of the energy-weighted sum-rule values and that 
the nucleus strongly absorbs energy from an external field, acting as a whole. Indeed, the 
giant resonances can be qualitatively described in terms of semiclassical hydrodynamical 
^ [ models PQ. 

Among many kinds of giant resonances, the isovector giant dipole resonance (GDR) is 
the most famous and exhausts almost 100 % of the sum-rule value. It is well explained 
by ordinary hydrodynamical models. This is rather exceptional because quantitative 
description of other modes requires a treatment as a Fermi liquid pQj. There are two 
famous hydrodynamical models for GDR: the Goldhaber- Teller (GT) model |2J and the 
Steinwedel- Jensen (SJ) model The GT model predicts A^ 1 ^ dependence of the GDR 
frequencies which arises from the concept that the restoring force is proportional to the 
nuclear surface area. In contrast, the SJ model relaxes the assumption of incompressibility, 
which leads to A" 1 / 3 dependence of the GDR frequencies. Experimental data are best 
fitted by a combination of these two [3]: In light nuclei, the data seem to indicate the A~ x ^ 
law, while the v4~ 1//3 dependence becomes increasingly dominant for increasing values of 
A. 

The hydrodynamical models have a close connection to the time- dependent Hartree- 
Fock (TDHF) theory. In the limit of h — > 0, the TDHF equation goes over into the Vlasov 
equation. Therefore, without the collision term, the TDHF should provide a microscopic 
description of an appropriate hydrodynamical model. Recently, we have proposed the 
real-time TDHF method combined with the absorbing-boundary condition (TDHF+ABC 
method) for a linear response function jolOj . This may be regarded as an extension of the 
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continuum random-phase approximation (RPA), made applicable to a deformed system. 
In this paper, we show real-time dynamics of the TDHF for GDR and discuss how their 
properties change from stable (N ~ Z) to unstable nuclei (N ^> Z). As is discussed 
in the following sections, we take a small-amplitude limit of the TDHF. Although this 
is equivalent to the RPA, the time-dependent snap shots of the TDHF wave packet may 
provide an intuitive dynamical picture of the GDR. It should be also noted that the TDHF 
provides a proper description for low-lying modes as well, for which quantum effects are 
so important that the semiclassical hydrodynamical models are not applicable. 

2. Linearized TDHF in real time 

The HF ground state is assumed to be a Slater determinant which consists of A single- 
particle orbitals, $q(xi, ■ ■ ■ ,Xa) = det{<f>i(xj)}i,j=i,-,A with x = (r,a,r). Each single- 
particle orbital is determined by 

h[(j), 4>*\<j)i{x) = ei4>i{x) for i = 1, • • • , A, (1) 

where /i[</>, </>*] is the single-particle Hamiltonian which depends on <pi{x) (i = 1, • ■ ■ ,A) 
self-consistently. The TDHF equation is obtained by replacing, in Eq. e,i by the 
time derivative ihd/dt, and <fii(x) by the time- dependent wave function ipi(x,t). The 
TDHF equation is now linearized with respect to variation of each single-particle wave 
function and a time-dependent external field v(x, t). Substituting i[)i(x,t) = (4>i(x) + 
5ipi(x,t))e~ leit / n into the TDHF equation, we have 

d 

ih—5il>i(x, t) = (h[<f>, 0*] - a) 5iJi(x, t) + 8h{t)^i{x) + v(x, t)<j>i{x), (2) 

where 8h(t) = h[if),if)*] — h[(p,4>*} in the first order of Sipi(x,t). If we put 8h(t) = 0, 
Eq. (j2J) gives unperturbed particle-hole excitations with a fixed single-particle potential 
in h[(j), (jf\. 5h(t) is nothing but the residual interaction in the language of the energy 
representation. The second term in the r.h.s. of Eq. (j2J) contains a dynamical effect 
which comes from variations of the self-consistent one-body potential. 

Equation (j2J) is equivalent to the well-known RPA equation in the energy representation. 
In practice, however, there are some differences, advantages and disadvantages in each 
method. For instance, the uncertainty in energy, AE, is inversely proportional to the 
period of the time propagation T; AE ~ h/T. Therefore, when we want to distinguish 
states nearly degenerate, we need to propagate the wave functions for a long period of 
time. In this case, the energy representation may be a better choice. On the other 
hand, when we are interested in a bulk structure of excited states in a wide range of 
energy, calculations using the time representation becomes more efficient than those with 
the energy representation. The time-dependent calculation should be suitable for giant 
resonances, since their energies are rather high and spread over a wide range of energy. 

The transition density in the time representation is defined by the density variation 
from its ground-state value, 

Sp(x,x';t) = p(x,x';t) — po(x,x') 

= jz{<PKx)5^{x\t) + 5^{x,t)U^)}- (3) 

8=1 
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In this paper, we are mainly interested in the spin- independent diagonal part of Eq. (J3J); 

A 

6p T {r; t) = E E {<£* t) + 6tP*(x, t)&(z)} . (4) 

1 = 1 c 

The expectation value of an operator F(f, r) can be expressed as 

F(t) = F + 6F(t) = F + E / d 3 rF(r, r)5p T (f; t), (5) 

where F is the ground-state expectation value. 

The external field v(x,t) in Eq. (J2J) can be chosen according to the purpose of the 
calculation. In order to calculate the strength distribution in a wide range of energy, an 
instantaneous external field, v(x, t) = v(x)S(t), is suitable, because this excites the system 
to states in all energies. In contrast, if we adopt an oscillating field with a fixed frequency 
u, v(x,t) = v(x) cos(ujt), the system is excited to a specific state with E x = fku. In this 
way, we can investigate dynamical properties of the specific state in the time-dependent 
manner. 

To calculate the strength function of the operator F(r, r), 
S(F; E) = E<H£ - E') |(^|F|$ )| 2 , (6) 

E' 

we adopt the external field v (x) proportional to F(r, r). Then, S(F; E) can be obtained 
as the Fourier transform of the expectation value of Eq. (j3J). 

Before showing results, let us discuss a numerical difficulty related to presence of zero 
modes in nuclei. The "zero mode" means zero-energy modes of excitation associated with 
the spontaneous symmetry breaking in the HF states, such as translation and rotation. 
These modes should correspond exactly to the zero energy if the numerical calculation is 
perfect. However, a small numerical error and approximation may give imaginary energies 
to these modes. Since the time evolution of wave functions carries all the information 
of the excited states, the presence of these imaginary-energy modes leads to a kind of 
numerical instability to prevent performing a long period of the time propagation. Thus, 
we need to remove components of the zero modes from the time- dependent single-particle 
wave functions 8ipi(x,t). The zero modes can be constructed by operating the symmetry 
operator P and its conjugate one Q to the ground state. For the translational case, P is 
the total momentum operator and Q is the center-of-mass coordinate. 



3. Giant dipole resonances in stable and neutron-rich nuclei 

We now apply the method to GDR in even-even Be isotopes. We calculated B(E1) 
distribution for Be isotopes in Ref. [H] using the SIII parameter set. We found that the 
large deformation splitting in 8 Be and 14 Be, because of the large quadrupole (prolate) 
deformation in the ground state (fy ~ 0.8). However, the width in 14 Be is much larger 
than that in 8 Be. These results are robust and do not depend on the choice of the Skyrme 
parameter set. We also found that there is a significant low-energy El strength around 
E x = 5 MeV for 14 Be However, the El strength and its peak position are rather 
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Figure 1. Calculated El oscillator strength distribution in ' Be. Thin solid and dashed 
lines indicate the response to dipole fields parallel and perpendicular to the symmetry 
axis, respectively. Thick line shows the total strength. 



sensitive to the choice of the parameters. Thus, in this section, we show time-dependent 
transition densities for the main peak of GDR. 

The calculation is performed on the three-dimensional Cartesian coordinate grid space, 
using the Skyrme energy functional of Ref. [7j. The Galilean symmetry is respected in 
this functional including the spin-orbit, Coulomb, and time-odd densities. We adopt the 
SGII parameter set in the calculation. In order to take account of the single-particle 
continuum, we use the absorbing potential outside of the interacting region 5 G 1 . In 
Fig.n we show the El oscillator strength distribution for 8 ' 14 Be. Two-peak structure due 
to the deformation splitting is prominent for both 8 Be and 14 Be. Hereafter, let us focus 
our discussion on the peak around E x = 15 MeV with K = 0. 

We use a Gaussian-pulse external field, v(x,t) = M(E1) K=0 cos(cjt)e _7 ^~* ^ , to selec- 
tively excite the GDR around E x = hco = 15 MeV, with 7 = 3 MeV/h and t = 2 h/MeV. 
Then, the spin-independent transition density of Eq. (@J) is calculated in the 3D coordinate 
space. It turns out that one of the Steinwedel- Jensen's assumptions, Sp n (r; t) = —5p p (f; t), 
is approximately satisfied for 8 Be. In contrast, in 14 Be, we see a large deviation from this 
property. Figure 121 shows how Sp T (f, t) (r = p, n) evolve in time in the x-z plane. The time 
difference from one panel to the next is At = 0.2 h/MeV which roughly corresponds to the 
half period tc/uj. We see that significant portions of neutrons actually move together with 
protons. The neutron transition density 5p n shows a peculiar node structure. In Fig. |21 
the regions of 5p p > (5p p < 0) have a large overlap with the those of 5p n > (5p n < 0). 
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Figure 2. Snap shots of calculated 6p T (r,t) in the x-z plane for the K = peak at 
i? x = 15 MeV in 14 Be. The upper panels (pl-4) indicate 5p p (f]t), while the lower (nl- 
4) for Sp n (f;t). White (black) regions indicate those of 5p T > (5p T < 0). The time 
difference between two neighboring panels is At = 0.2 h/MeV. The two panels at the 
same column corresponds to the same time t. 

This means a violation of the property of the SJ model, Sp p + 5p n = 0. Weakly bound 
neutrons in neutron-rich nuclei seem to be significantly affected by strong attraction be- 
tween protons and neutrons, and to oscillate in phase with protons' movement. This is a 
consequence of the dynamical effect of the time-dependent self-consistent potential, 5h(t) 
inEq. ©. 

4. Summary 

We present calculations of the linearized TDHF method in real time for giant dipole 
resonances in 8 ' 14 Be. These nuclei are calculated to be largely deformed in the ground 
state. The main GDR peak is split into two peaks with K = and K = 1. We show 
the time-dependent transition densities for K = peaks. The total density, p p + p n , is 
approximately conserved for 8 Be, while its conservation is significantly violated in 14 Be. 
The time evolution of the transition density, 5p T (t), suggests a strong dynamical effect for 
neutron-rich nuclei, and seems to indicate the mixture of the isoscalar volume- type and 
the isovector surface-type components. 
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6 



REFERENCES 

1. P. Ring and P. Schuck, The Nuclear Many-Body Problem ( Springer- Verlag, New York, 
1980). 

2. M. Goldhaber and E. Teller, Phys. Rev. 74 (1948) 1046. 

3. H. Steinwedel and J. H. D. Jensen, Z. Naturforschung 5A (1950) 413. 

4. B. L. Berman and S. C. Fultz, Rev. Mod. Phys. 47 (1975) 713. 

5. T. Nakatsukasa and K. Yabana, J. Chem. Phys. 114 (2001) 2550. 

6. T. Nakatsukasa and K. Yabana, Phys. Rev. C 71 (2005) 024301. 

7. P. Bonche and H. Flocard and P. H. Heenen, Nucl. Phys. A 467 (1987) 115. 



